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THEORY  OF  IMPURITY-SHIFTED  INTERSUBBAND  TRANSITIONS  IN 

N-TYPE  INVERSION  LAYERS  ON  (100)  SILICON* 

G.  M.  Kramer  and  R.  F.  Wallis 

University  of  California,  Irvine,  California  92717,  USA 

Binding  energies  and  variational  wave  functions  are 
calculated  as  functions  of  electric  field  for  elec- 
tron impurity  states  split  off  from  electric  subbands 
of  (100)  silicon  and  silicon  dioxide.  The  impurity- 
shifted  intersubband  transition  energies  are  pre- 
sented. 

1.  Introduction 

The  development  of  MOS  devices  has  stimulated  a great  deal 
of  interest  in  the  basic  properties  of  inversion  layers  at 
semiconductor  boundaries.  Of  some  importance  are  impuri- 
ties and  their  associated  electronic  bound  states.  We 
have  made  a theoretical  investigation  of  several  bound 
states  split  off  from  the  lowest  three  electric  subbands 
of  an  n-type  inversion  layer  on  (100)  silicon  bounded  by 
silicon  dioxide. 


2.  Theoretical  Development 

We  consider  two  contiguous,  semi-infinite  half-spaces,  one 
of  p-type  silicon  and  the  other  of  silicon  dioxide,  with 
a common  boundary  parallel  to  a (100)  plane  of  the  silicon. 
We  assume  that  the  potential  energy  of  an  electron  under- 
goes an  infinite  discontinuous  jump  as  the  electron  passes 
from  the  silicon  into  the  silicon  dioxide.  In  the  region 
of  the  inversion  layer,  there  is  an  electric  field  which, 
to  a good  first  approximation?1  can  be  taken  to  be  constant. 

We  assume  that'  an  impurity  ion  of  charge  +Ze  Is  located  at 

the  boundary  between  the  silicon  and  the  oxide.  For  this 
2 

situation  , the  lowest  lying  bound  states  associated  with 
a (100)  Interface  on  silicon  are  derived  from  the  constant 
energy  ellipsoids  whose  major  axes  are  perpendicular  to  the 
interface. 


The  Hamiltonian  is  taken  to  be  
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H.  • - —^5  - — — * + U(r)  , (3) 

1 ax2  ay2 

Y - 5 " (€2“€1)/4Z€2»  H-fi4(614^2)3/4Z3mZe5,  C is 

the  electric  field  in  esu,  mt  and  m^  are  the  transverse  and 
longitudinal  effective  masses,  and  and  €2  are  the  diel- 
ectric constants  of  S102  and  silicon,  respectively.  In 
Eq.  (1),  energy  is  measured  in  units  of  the  effective 
Rydberg t Ry*  - 2mte4Z2/h2(€1+€2)2,  and  length  in  units  of 
the  effective  Bohr  radius,  a*  - h2(61+^0)/2Zm.e2.  For  the 
slllcon-S102  system,  Ry  >42.3  meV  and  a*  - 21.8  A.  The 
impurity  ion  is  taken  to  be  at  the  sillcon-S102  Interface 
and  screening  is  neglected,  so  U(r)  has  the  Coulomb  form 


U(r).--§-  , (4) 

M 

The  term  involving  6 in  Eq.  (2)  represents  the  interaction 
of  the  electron  with  its  own  image. 


In  the  limit  of  high  electric  fields,  an  adiabatic  theorem 
applies  and  the  impurity  wave  function  has  the  form 

♦<r)  - X^n) (x»y)*n<z)  (5) 

where  fQ(z)  is  an  eigenfunction  of  HQ  given  by  Eq.  (2)  and 
corresponds  to  the  electric  subband  of  interest.  We  con- 
sider three  subbands  and  choose  the  f (z)  to  be  variation-  ■ 
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al  functions  of  the  Fang-Howard  type 

fQ(z)  - Aqz  exp(-JbQz)  (6a) 

fj^z)  - AjZ(l -orjZjexpC-ibjZ)  (6b) 

f2(z)  - A2z(1  +or2z  +e2z2)exp(-ib2z)  . (6c) 

The  quantities  bQ  are  variational  parameters,  the  Aq  are 
normalization  constants,  and  or^,  cr2,  and  02  are  chosen  so 
that  the  functions  fQ(z)  are  orthogonal  on  the  Interval 
0<z<«,  We  take  to  have  either  jaf  _two  forma  

Xon)(*»7)  - BQ  exp[-$a<n)(x2 +y2)*]  (7a) 


(7b) 


X^U.y)  - B^xiiyjwcpC-i  a<n)(x2  + y2)*] 

(n\  3 

where  the  B,  are  normalization  constants  and  the  a;  ' are 
variational  parameters.  The  parameters  a£  ' and  bQ  are 
varied  independently  to  minimize  the  expectation  value  of 
the  Hamiltonian,  <H>.  The  expectation  value  of  HQ,  < HQ> , 
is  also  minimized  separately  with  respect  to  the  bQ  to  give 
the  continuum  energy.  The  impurity  binding  energy,  Eg,  is 
then  given  by 

Eg  - <H0>  - <H>  . (8) 

3.  Results  and  Discussion 

The  Impurity  binding  energies  have  been  calculated  for 
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electric  fields  between  10  and  10  esu  for  the  three  low- 
est subbands  and  m-0,  ±1.  The  results  are  plotted  in  Fig.  1. 


Fig.  1.  Impurity  state  binding  energies  versus  field 

The  Impurity  binding  energy  decreases  as  one  goes  from  the 
lowest  subband  to  successively  higher  subbands  for  both 
the  m-0  and  m-*l  bound  states.  For  a given  bound  state 
associated  with  a given  subband,  the  binding  energy  in- 
creases with  increasing  electric  field  and  approaches  the 
two-dimensional  value  in  the  limit  of  infinite  field. 

Electric  dipole  transitions  between  subbands  are-allowed- 
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transitions.  In  Fig.  2 are  given  our  calculated  impurity 
shifted  transition  energies  for  the  0-1  and  0-2  inter- 


subband  transitions  when.. 
o-O . The  transition  ener- 
gies are  shifted  to  higher 
energies  when  the  electron 
becomes  localized  at  the 

impurity.  At  a field  of 
2 

10  esu,  the  impurity 
shifted  0-1  transition  is 
70  percent  higher  than  the 
unshlfted  transition  and 
is  only  15  percent  lower 
than  the  0-2  unshlfted 
transition.  The  latter  re- 
sult is  in  qualitative 
agreement  with  recent  ex- 
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perimental  data.  As  the 
electric  field  increases, 
the  impurity  shifted  0-1 
transition  drops  back 
toward  the  0-1  free  elec- 
tron transition  and  away  fr 
sit ion. 


and  free  states 

the  0-2  free  electron  tran- 
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